Inhomogeneous Superconductivity in Comb-Shaped Josephson Junction Networks 
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We show that some of the Josephson couplings of junctions arranged to form an inhomogeneous 
network undergo a non-perturbative renormalization provided that the network's connectivity is 
pertinently chosen. As a result, the zero-voltage Josephson critical currents I c turn out to be 
enhanced along directions selected by the network's topology. This renormalization effect is possible 
only on graphs whose adjacency matrix admits an hidden spectrum (i.e. a set of localized states 
disappearing in the thermodynamic limit). We provide a theoretical and experimental study of this 
effect by comparing the superconducting behavior of a comb-shaped Josephson junction network 
and a linear chain made with the same junctions: we show that the Josephson critical currents 
of the junctions located on the comb's backbone are bigger than the ones of the junctions located 
on the chain. Our theoretical analysis, based on a discrete version of the Bogoliubov-de Gennes 
equation, leads to results which are in good quantitative agreement with experimental results. 
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It is a common belief that Josephson Junction Net- 
works (JJN) may be regarded as the prototype of a com- 
plex physical system with a variety of interesting physical 
behaviors, adjustable acting only on a few external pa- 
rameters and, by means of the modern fabrication tech- 
nologies, also on the building topology and geometry of 
the array [1]. It is now possible to build arrays with 
very small junctions [2] to detect effects due to single 
electrons in a range of temperatures related to another 
relevant energy scale: the charging energy. The possi- 
bility of controlling experimentally the competition be- 
tween Josephson and charging energies makes JJNs use- 
ful devices to investigate quantum phase transitions [3], 
or to model the physical properties of some real materi- 
als like granular superconductors [4] . Many of the results 
valid for JJNs are shared by cold atoms in optical lattices 
[5] since, in these systems, bosonic Josephson junctions 
and arrays may be rather easily realized [6]; in addition, 
Josephson networks and devices pave a very promising 
avenue to the quantum engineering of states relevant for 
quantum computing [7]. 

Inhomogeneous superconducting networks have been 
studied [8] mainly to provide a better understanding of 
the properties of well controlled disordered granular su- 
perconductors [9]. The appealing perspective to realize 
devices for the manipulation of quantum information re- 
cently stimulated the analysis of inhomogeneous planar 
JJNs with non conventional connectivity [10], engineered 
to sustain a topologically ordered ground state [11]. 
Transport measurements on superconducting wire net- 
works evidenced - in a pure system with non-dispersive 
eigenstates- interesting anomalies of the network critical 



current induced by the interplay between the network's 
geometry and topology and an externally applied mag- 
netic field [12]; more recently, the theoretical analysis of 
rhombi chains has evidenced the exciting possibility of 
being able to detect 4e superconductivity through mea- 
surements of the supercurrent in presence of a pertinent 
external magnetic field [13]. 

In this Letter we show that, even in absence of an 
externally applied magnetic field, a JJN fabricated on a 
pertinent graph [14] may support anomalous behaviors of 
the Josephson critical currents, which are induced by a 
non-perturbative renormalization of some of the Joseph- 
son couplings of the array. Our analysis clearly evidences 
that this renormalization is only attainable for the class 
of graphs, whose adjacency matrix supports an hidden 
spectrum [14,15]; thus, our findings are not generic to any 
inhomogeneous JJN. For instance, in absence of an exter- 
nal magnetic field, the networks analyzed in [8-10,12,13] 
should not give rise to any of the anomalous behaviors of 
I c discussed in this paper. 

In the following, we provide a theoretical and exper- 
imental study of the behavior of the Josephson critical 
currents measurable in a comb-shaped JJN made of Nb 
grains located at the vertices of a "comb" graph and 
linked by Josephson junctions [see Fig.l]. We compare 
our results with those obtained for a linear Josephson 
junction chain fabricated with the same junctions. Since 
one may regard the backbone of a comb graph as a dec- 
orated chain, it appears natural to compare its super- 
conducting properties with those of a linear chain since 
the latter is the simplest network with euclidean dimen- 
sion one. The result of this comparison shows that the 
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Josephson critical currents of the junctions located on 
the comb's backbone are sensibly bigger than the ones of 
the junctions located on the chain. 

Another way to look at a comb-shaped JJN is to re- 
gard it as a linear chain immersed in an environment 
mimicked by the addition of the fingers [16]. As in many 
Josephson devices one should then expect that the nom- 
inal value of the Josephson energy Ej of the junctions in 
the array gets renormalized by the interaction with the 
environment. This situation is often analyzed using ei- 
ther the Caldeira-Leggett [17] or the electromagnetic en- 
vironment [18] models. In these approaches one usually 
assumes that the effective boundary conditions for the 
quantum fluctuations of the environment modes do not 
depend on the Josephson couplings or on the network's 
topology: while this assumption is perfectly legitimate 
for weak environmental fluctuations, better care should 
be used if these fluctuations are strong as it may well 
happen for one dimensional JJNs. A simple paradigmatic 
example of a non perturbative renormalization of Joseph- 
son couplings is given by the simple inhomogeneous one- 
dimensional array analyzed in [19,20], where the source 
of inhomogeneity is given by putting on a site of the lin- 
ear chain a test junction with a different nominal value 
of the Josephson coupling Ej. In the sequel we show 
that, for a comb-shaped JJN, the Josephson couplings 
on the backbone get renormalized. Our explicit com- 
putation shows that this renormalization is indeed non 
perturbative since the peculiar connectivity of a comb 
modifies the spectrum of quantum modes living on lin- 
ear chains by the (obviously non-perturbative) addition 
of an infinite set of localized states, which disappear in 
the thermodynamic limit (the hidden spectrum): adding 
the fingers to a backbone chain is, in fact, a topolog- 
ical operation since it amounts to a non trivial change 
of boundary conditions for the Josephson linear chain. 
In a different context, the interplay between an hidden 
spectrum and a change in boundary conditions has been 
recently used in [21]. 

We use the lattice Bogoliubov-de Gennes (LBdG) 
equations [22] to compare the properties exhibited by 
Josephson linear chains and comb-shaped Josephson net- 
works. Using the eigenfunctions of the LBdG equations, 
a self-consistent computation yields for both systems the 
gap function, the chemical potential and the quasi par- 
ticle spectrum. We show that, for a linear chain, the 
superconducting gap and critical temperature satisfy to 
the well-known BCS equations and that, on the backbone 
of a comb JJN, the BCS equations are satisfied with a 
renormalized value of the Josephson energy. Then, we 
compute the zero-voltage Josephson critical currents I c 
on the comb's backbone and compare our results for I c 
with the outcomes of experimental measurements: our 
computation not only confirms with good accuracy the 
experimental results of [23], but is also in good agree- 
ment with new data obtained at temperatures closer to 
the critical temperature for the onset of superconductiv- 
ity in Nb grains. The new data are shown in Fig. 2. 



To obtain a discrete version of the BdG equations suit- 
able to describe the JJNs fabricated in [23] , we make the 
ansatz that the eigenfunctions of the continuous BdG 
equations [22] may be written in a tight binding form as 
u a (r) = Y,iU a (i)4>i(r) and v a (r) = J^i v a (i)4>i(r); i la- 
bels the position of a superconducting island while the 
contribution of the electronic states participating to su- 
perconductivity in a given island is effectively described 
by a field 4>i(r), whose specific form depends on the ge- 
ometry of the islands and on the fabrication parameters 
of the connecting junctions. The LBdG equations then 
read 

€ Q u a (i) = y^e»jMa(j) + A(i)v a (i) (1) 
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G a v a (i) = - ejjV a (j) + A*(i)u a (i). (2) 

3 

where u a and v a satisfy to J2i [l u a(*)| 2 + kaWI 2 ] = 1- 
The matrix is defined by = — tA^j + U(i)8ij — 
fiSij, with Aij being the adjacency matrix of the net- 
work [24], fx = [i - J df(j)i(r} (-% 2 V 2 /2m) <f>i(r) and 
t = - J df(j3 2 (r)[~h 2 V 2 /2m + C/ (r)]^(r) » Ej. Ej = 
(h/2e)I c is the nominal value of the Josephson energy 
of all the junctions in the network while I c is the un- 
renormalized zero-voltage Josephson critical current of 
each junction. Uo(r) mimics the effects of the barrier 
between the superconducting islands. Self-consistency 

requires A(i) = VEa iiaB^W^If 6«J an d 
U{i) = -V£«[KW| 2 /a + K(i)| 2 (1 -/«)], where 
V = V(f) 2 (r = fi) is assumed to be independent on i. 
Topology is encoded in the term —tAij appearing in the 
definition of the matrix e^ , while the specific values of 
t and V depend - as a result of our ansatz on the form 
of the eigenfunctions of the BdG equations- only on the 
Mr)- 
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FIG. 1. Schematic drawing of a comb array. The super- 
conducting islands (full box) are connected in series to each 
other through Josephson junctions (JJs). The finger arrays 
are connected to each other only through JJs to the central 
islands forming the backbone array. 

To justify the assumptions involved in the derivation 
of Eqs.(l)- (2), we observe that, for the JJNs described 
in [23], capacitive (inter islands and with a ground) ef- 
fects are negligible, that the total number of electrons on 
the island Af is much larger than the number of electrons 
tunneling through the Josephson junction and that all 
islands contain approximately the same Af (Af(i) = Af). 
Furthermore, the islands are big enough to support the 
same superconducting gap of the Nb bulk material. As a 
result one may require (j>i(f) to be position-independent 
on each island except for a small region near the junc- 
tion and to be the same on each island with a nor- 
malization given by J dr<f>i(r)<f>i{f) = Af(i) = Af and 
/ dr(f>i(f)(j)j(r) w for i ^ j. In our derivations we put 
Af =1. 

For a linear array, the LBdG may be readily solved 
leading to an uniform potential U(i) = U c and an uniform 
pair potential A(z) = A c . From the eigenvalue equation 
-EjY.jAijifkii) = efeV'fe(i), one gets e fe = —2Ejcosk: 
it follows = a/ A 2 + Ef. with E k = e k + U c - fi. The 
BCS-like behaviour is obtained when Ek = ek, which 
happens since U c — and fi = E F . When A c /Ej <C I, 

for T = 0, one gets A C (T = 0) = 8Eje~ 2nEj ^, while, 
for T = T C (i.e., A C (T = T c ) = 0), one obtains k B T c = 
CEje- 2nEj ^ ', with C = 4.54. It is comforting that the 
assumptions on which our approach is based lead, for 
the chain, to results having the same functional form of 
the well-known BCS formulas for the gap at T — (i.e., 
A(T = 0) = 2huj D e- 1 / n ^ VBCS ) and the BCS critical 
temperature (i.e., k B T c = l.UhuDe' 1 ^^^ 03 ), pro- 
vided that ti(0)Vbcs "C 1 [22]: in addition, one gets also 
A C (T = 0)/k B T c = 8/C w 1.76. 

Measurements on a chain made with Nb grains yield 
T c w 8.8K and A C (T = 0) w lAmeV w k B ■ 15. 9 K; 
furthermore, in the experimental setup described in [23] 
it is I c 18/iA. The parameters Ej and V, deter- 
mined from the BCS equation yielding the chain's criti- 
cal temperature, are then given by Ej w k B ■ 430K and 
V/Ej = 1.185. In Fig. 2 we plot for several temperatures 
the measured I c (circles) and the critical currents ob- 
tained inserting A C (T) in the well-known Ambegaokar- 
Baratoff expression [25] for the zero-voltage Josephson 
current (lower solid curve): the agreement is excellent. 

For a comb network with N x N islands (see Fig. I), 
one finds a solution of the LBdG equations (I)- (2) 
where both the Hartree-Fock potential U(i) and the 
gap function A(i) are position dependent. We denote 
the islands by (x,y), x labeling the finger and \y\ the 
distance from the backbone, expressed in lattice units. 
The eigenvalue equation —EjJ^j Aijip a (j) = e a ip a (i), 
admits [24], in addition to a set of delocalized states 



with energies ranging from —2Ej to 2Ej, a localized 
ground-state tp = (C /V~N)e~ y ^ , corresponding to the 
eigenvalue e = —2\/2Ej (C 2 = 1/V2 and £ given by 
sinh(l/£) = 1) and an hidden spectrum made of other 
eigenstates localized around the backbone [24]. For a 
crude analytical estimate, one may require that, away 
from the backbone, the fingers may be regarded as a 
linear chain with uniform potentials (i.e., A(i) = A c 
and U(i) — U c ). To get then coupled equations for A&, 
A c , Ub, and U c , one writes the LBdG equations (l)-(2) 
on a backbone's grain i. We set u a (i) = U a ip a (i) and 
v a (i) — V a tf a (i), with U 2 + V 2 = 1. The self-consistency 

vc 2 

equation for U implies that, at T = 0, Ub ~ U c 

upon requiring fx w Ub one immediately sees that, due to 
the localized modes in the fcrmionic spectrum, the chem- 
ical potential on the comb's backbone is smaller than the 
one measured on the chain. 

By substituting the wavefunctions of the eigenstates of 
the hidden spectrum [24] in Eqs.(l)-(2) and using jl w Ub 
one gets 

A b V r /2 j; cosfc //3 \ 

A h = A c + / dk • tanh £ e k . 

n Jo e fe VI + cos 2 k V 2 / 

(3) 

where e fe = ^A 2 + AEj (1 + cos 2 k). The hidden spec- 
trum eigenstates contribute to the gap function A^ 
through the second term in the rhs of Eq.(3): without 
them, A;, equals A c . 

When Ej > A 6 ,A C , Eq.(3), at T = 0, yields 

OTuT = 7^ = K where 1° = log (1 + V2) ■ 

2*Ej 

Furthermore, at low temperatures, Af,(T)/A c (T) w 
A 6 (T = 0)/A c (T = 0). Using the parameters Ej and V 
obtained from the measurements carried on the JJ chain, 
for a JJ comb one gets JC w 1.13. 

Upon requiring that, as for the linear chain, the T = 
backbone's gap function has a BCS like functional form, 

i.e. A 6 (T = 0) = 8E J e- 2 ' rEj ^, with Ej and V the 
renormalized Josephson energy and the renormalized in- 
teraction term, one is able to estimate the renormal- 
ization of the Josephson coupling within the LBdG ap- 
proach. Namely, one has, 

Ej = JCEj; V = KV, (4) 

which embodies the effects of the hidden spectrum on the 
Josephson couplings. 

In Fig. 2 we plot, as a function of the normalized tem- 
perature, the values of I c measured with the methods 
described in [23] (squares) and the values of I c obtained 
from the Ambegaokar-Baratoff formula using both the 
renormalized coupling given by Eq.(4) and the gap func- 
tion along the backbone for the comb-like JJN studied 
in [23] (solid curve): the agreement between theory and 
experiments is very good at low temperatures, while the 
theory gives a slight overestimate at higher temperature. 
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FIG. 2. Critical currents (in units of the critical current on 
the reference chain at T = 1.2K) as a function of T/T c for 
the backbone and the chain. The solid lines are the estimated 
critical currents for the backbone (top) and the chain (bot- 
tom). Circles (squares): experimental values for the chain 
(backbone). 

We showed that a non perturbative (i.e. induced by the 
states of the hidden spectrum) renormalization of some of 
the Josephson couplings of a comb-shaped JJN is respon- 
sible for the observed enhancement of I c of the Josephson 
junctions located along the comb's backbone. The key 
assumption in our derivation is that the eigenfunctions 
of the BdG equations may be written in a tight binding 
form; once this assumption is made, one is able to derive 
Eqs. (1)- (2) and to account for all the dependence on 
the electronic states into the definition of the parameters 
Ej and V, which, in this paper, we determined from the 
measurements carried on the linear chain. Our approach 
yields a value of the renormalized Josephson coupling of 
the junctions located on the comb's backbone in excel- 
lent agreement with the experimental results (see fig. 2). 
We expect that similar phenomena happen for the class 
[15] of JJNs fabricated on graphs whose adjacency matrix 
supports an hidden spectrum. 
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